Comparing surface-soil moisture from the SMOS mission and the ORCHIDEE land-surface model over the Iberian Peninsula

The aim of this study is to compare the surface soil moisture (SSM) retrieved from ESA's Soil Moisture and Ocean Salinity mission (SMOS) with the output of the ORCHIDEE (ORganising Carbon and Hydrology In Dynamic EcosystEm) land surface model forced with two distinct atmospheric data sets for the period 2010 to 2012. The comparison methodology is first established over the REMEDHUS (Red de Estaciones de MEDición de la Humedad def Suelo) soil moisture measurement network, a 30 by 40. km catchment located in the central part of the Duero basin, then extended to the whole Iberian Peninsula (IP). The temporal correlation between the in-situ, remotely sensed and modelled SSM are satisfactory (r. >. 0.8). The correlation between remotely sensed and modelled SSM also holds when computed over the IP. Still, by using spectral analysis techniques, important disagreements in the effective inertia of the corresponding moisture reservoir are found. This is reflected in the spatial correlation over the IP between SMOS and ORCHIDEE SSM estimates, which is poor (¿. ~. 0.3). A single value decomposition (SVD) analysis of rainfall and SSM shows that the co-varying patterns of these variables are in reasonable agreement between both products. Moreover the first three SVD soil moisture patterns explain over 80% of the SSM variance simulated by the model while the explained fraction is only 52% of the remotely sensed values. These results suggest that the rainfall-driven soil moisture variability may not account for the poor spatial correlation between SMOS and ORCHIDEE products.Peer ReviewedPostprint (published version

Crop Yield Estimation and Interpretability With Gaussian Processes

This work introduces the use of Gaussian processes (GPs) for the estimation and understanding of crop development and yield using multisensor satellite observations and meteo- rological data. The proposed methodology combines synergistic information on canopy greenness, biomass, soil, and plant water content from optical and microwave sensors with the atmospheric variables typically measured at meteorological stations. A com- posite covariance is used in the GP model to account for varying scales, nonstationary, and nonlinear processes. The GP model reports noticeable gains in terms of accuracy with respect to other machine learning approaches for the estimation of corn, wheat, and soybean yields consistently for four years of data across continental U.S. (CONUS). Sparse GPs allow obtaining fast and compact solutions up to a limit, where heavy sparsity compromises the credibility of confidence intervals. We further study the GP interpretability by sensitivity analysis, which reveals that remote sensing parameters accounting for soil moisture and greenness mainly drive the model predictions. GPs finally allow us to identify climate extremes and anomalies impacting crop productivity and their associated drivers

Impact of day/night time land surface temperature in soil moisture disaggregation algorithms

Since its launch in 2009, the ESA’s SMOS mission is providing global soil moisture (SM) maps at ~40 km, using the first L-band microwave radiometer on space. Its spatial resolution meets the needs of global applications, but prevents the use of the data in regional or local applications, which require higher spatial resolutions (~1-10 km). SM disaggregation algorithms based generally on the land surface temperature (LST) and vegetation indices have been developed to bridge this gap. This study analyzes the SM-LST relationship at a variety of LST acquisition times and its influence on SM disaggregation algorithms. Two years of in situ and satellite data over the central part of the river Duero basin and the Iberian Peninsula are used. In situ results show a strong anticorrelation of SM to daily maximum LST (R˜-0.5 to -0.8). This is confirmed with SMOS SM and MODIS LST Terra/Aqua at day time-overpasses (R˜-0.4 to -0.7). Better statistics are obtained when using MODIS LST day (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.06 m3 /m3 ) than LST night (R˜0.45 to 0.80; ubRMSD˜0.04 to 0.07 m3 /m3 ) in the SM disaggregation. An averaged ensemble of day and night MODIS LST Terra/Aqua disaggregated SM estimates also leads to robust statistics (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.07 m3 /m3 ) with a coverage improvement of ~10-20 %.Peer ReviewedPostprint (published version

L-Band Vegetation optical depth and effective scattering albedo estimation from SMAP

Over land the vegetation canopy affects the microwave brightness temperature by emission, scattering and attenuation of surface soil emission. Attenuation, as represented by vegetation optical depth (VOD), is a potentially useful ecological indicator. The NASA Soil Moisture Active Passive (SMAP) mission carries significant potential for VOD estimates because of its radio frequency interference mitigation efforts and because the L-band signal penetrates deeper into the vegetation canopy than the higher frequency bands used for many previous VOD retrievals. In this study, we apply the multi-temporal dual-channel retrieval algorithm (MT-DCA) to derive global VOD, soil moisture, and effective scattering albedo estimates from SMAP Backus-Gilbert enhanced brightness temperatures posted on a 9 km grid and with three day revisit time. SMAP VOD values from the MT-DCA follow expected global distributions and are shown to be highly correlated with canopy height. They are also broadly similar in magnitude (though not always in seasonal amplitude) to European Space Agency Soil Moisture and Ocean Salinity (SMOS) VOD. The SMOS VOD values are based on angular brightness temperature information while the SMAP measurements are at a constant incidence angle, requiring an alternate approach to VOD retrieval presented in this study. Globally, albedo values tend to be high over regions with heterogeneous land cover types. The estimated effective scattering albedo values are generally higher than those used in previous soil moisture estimation algorithms and linked to biome classifications. MT-DCA retrievals of soil moisture show only small random differences with soil moisture retrievals from the Baseline SMAP algorithm, which uses a prior estimate of VOD based on land cover and optical data. However, significant biases exist between the two datasets. The soil moisture biases follow the pattern of differences between the MT-DCA retrieved and Baseline-assigned VOD values

Multi-temporal evaluation of soil moisture and land surface temperature dynamics using in situ and satellite observations

Soil moisture (SM) is an important component of the Earth’s surface water balance and by extension the energy balance, regulating the land surface temperature (LST) and evapotranspiration (ET). Nowadays, there are two missions dedicated to monitoring the Earth’s surface SM using L-band radiometers: ESA’s Soil Moisture and Ocean Salinity (SMOS) and NASA’s Soil Moisture Active Passive (SMAP). LST is remotely sensed using thermal infrared (TIR) sensors on-board satellites, such as NASA’s Terra/Aqua MODIS or ESA & EUMETSAT’s MSG SEVIRI. This study provides an assessment of SM and LST dynamics at daily and seasonal scales, using 4 years (2011–2014) of in situ and satellite observations over the central part of the river Duero basin in Spain. Specifically, the agreement of instantaneous SM with a variety of LST-derived parameters is analyzed to better understand the fundamental link of the SM–LST relationship through ET and thermal inertia. Ground-based SM and LST measurements from the REMEDHUS network are compared to SMOS SM and MODIS LST spaceborne observations. ET is obtained from the HidroMORE regional hydrological model. At the daily scale, a strong anticorrelation is observed between in situ SM and maximum LST (R ˜ -0.6 to -0.8), and between SMOS SM and MODIS LST Terra/Aqua day (R ˜ - 0.7). At the seasonal scale, results show a stronger anticorrelation in autumn, spring and summer (in situ R ˜ -0.5 to -0.7; satellite R ˜ -0.4 to -0.7) indicating SM–LST coupling, than in winter (in situ R ˜ +0.3; satellite R ˜ -0.3) indicating SM–LST decoupling. These different behaviors evidence changes from water-limited to energy-limited moisture flux across seasons, which are confirmed by the observed ET evolution. In water-limited periods, SM is extracted from the soil through ET until critical SM is reached. A method to estimate the soil critical SM is proposed. For REMEDHUS, the critical SM is estimated to be ~0.12 m3/m3 , stable over the study period and consistent between in situ and satellite observations. A better understanding of the SM–LST link could not only help improving the representation of LST in current hydrological and climate prediction models, but also refining SM retrieval or microwave-optical disaggregation algorithms, related to ET and vegetation status.Peer ReviewedPostprint (published version

Comparison of measured brightness temperatures from SMOS with modelled ones from ORCHIDEE and H-TESSEL over the Iberian Peninsula

L-band radiometry is considered to be one of the most suitable techniques to estimate surface soil moisture (SSM) by means of remote sensing. Brightness temperatures are key in this process, as they are the main input in the retrieval algorithm which yields SSM estimates. The work exposed compares brightness temperatures measured by the SMOS mission to two different sets of modelled ones, over the Iberian Peninsula from 2010 to 2012. The two modelled sets were estimated using a radiative transfer model and state variables from two land-surface models: (i) ORCHIDEE and (ii) H-TESSEL. The radiative transfer model used is the CMEM. Measured and modelled brightness temperatures show a good agreement in their temporal evolution, but their spatial structures are not consistent. An empirical orthogonal function analysis of the brightness temperature's error identifies a dominant structure over the south-west of the Iberian Peninsula which evolves during the year and is maximum in autumn and winter. Hypotheses concerning forcing-induced biases and assumptions made in the radiative transfer model are analysed to explain this inconsistency, but no candidate is found to be responsible for the weak spatial correlations at the moment. Further hypotheses are proposed and will be explored in a forthcoming paper. The analysis of spatial inconsistencies between modelled and measured TBs is important, as these can affect the estimation of geophysical variables and TB assimilation in operational models, as well as result in misleading validation studies